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Abstract 

We combine the six-dimensional helicity formalism of Cheung and O'Connell with Z?-dimensional 
generalized unitarity to obtain a new formalism for computing one-loop amplitudes in dimensionally 
regularized QCD. With this procedure, we simultaneously obtain the pieces that are constructible 
from four-dimensional unitarity cuts and the rational pieces that are missed by them, while retain- 
ing a helicity formalism. We illustrate the procedure using four- and five-point one-loop amplitudes 
in QCD, including examples with external fermions. We also demonstrate the technique's effective- 
ness in next-to-leading order QCD corrections to Higgs processes by computing the next-to-leading 
order correction to the Higgs plus three positive-helicity gluons amplitude in the large top-quark 
mass limit. 
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I. INTRODUCTION 



With the LHC recently surpassing the Tevatron for the world record in hadron collider 
beam intensity and energy, we stand ready to test and expand our current understanding of 
particle physics. The LHC is well on its way in hunting for the Higgs boson and for physics 
beyond the Standard Model. To engage in these ambitious searches, however, a detailed 
picture of Standard Model backgrounds is necessary. Higher order QCD is an important 
theoretical tool in both Higgs production and the Standard Model background. 

The QCD corrections are large and much effort has been dedicated to calculating im- 
portant cross sections [1]. While some processes will require next-to- next-to-leading order 
precision and resummation of large logarithms, for many processes knowledge through next- 
to-leading order (NLO) in the perturbative expansion will be sufficient. This need for NLO 
calculations has been discussed and codified in an experimenter's wishlist jlj. One-loop 
calculations in QCD are therefore a crucial part of theoretical studies of backgrounds and 
signals at the LHC. 

As the number of external partons increases, Feynman-diagrammatic techniques en- 
counter difficulties, especially at loop level. This has motivated many attempts to seek 
improved techniques. In recent years there has been rapid progress significantly reducing 
the computational intensity behind NLO calculations. The unitarity method, pioneered by 
Bern, Dixon, Dunbar and Kosower in the mid-nineties 0, [sj] , offers an alternative using the 
branch-cut structure of loop integrals to find the coefficients of the integrals in terms of 
products of on-shell tree amplitudes. Use of generalized unitarity is natural in this frame- 
work, giving further enhancements to the method [ij. The tree amplitudes are sewn together 
across various generalized channels to reconstruct the loop amplitude. In this way only phys- 
ical degrees of freedom are used. One-loop amplitudes for various supersymmetric theories 
can be constructed in their entirety through the use of four- dimensional cuts, though non- 
supersymmetric theories^including QCD, have cut-free rational terms that must be handled 
through other means 0, HJ. By using complex momenta, Britto, Cachazo and Feng showed 
that generalized quadruple cuts determine all scalar box integral coefficients as products 
of four on-shell tree amplitudes j6[, motivating new reduction procedures on the integrand. 
Ossola, Papadopoulos and Pittau (OPP) followed suit with a purely algebraic method to 
reduce tensor integrals to a basis of master scalar integrals 0] in a form directly usable in 
numerical programs. The OPP method uses a particular parametrization of the loop mo- 
mentum to find coefficients for all scalar integrals by solving a system of linear equations. 
Forde followed with a method for simple extraction of triangle and bubble coefficients in 
four- dimensional cuts that exploits the analytic structure of the integrand [&]. Coefficients 
for all scalar integrals can be found through the limiting behavior of the products of tree 
amplitudes by using the momentum parametrization in ref. Q, first given in ref. joj. A 
numerical variant of Forde's method was then developed and applied in ref. llOj. There have 
also been important improvements in Feynman diagrammatic techniques [111 ] . 

Another important technique that greatly simplifies amplitude calculations is helicity 
methods [l2[. At tree level these lead to enormous simplifications, especially for massless 
amplitudes. The massless spinor-helicity formalism in four dimensions has been an invaluable 
tool for simplifying calculations and exposing the structure of amplitudes starting with the 
discovery of the Parke- Taylor formula for maximally helicity- violating tree amplitudes (l3| . 
At loop level helicity methods have also played an important role in pushing back the 
calculational frontiers. Recently the first NLO calculation of a hadron-collider process with 
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five final-state objects was presented: W+4-jet production [14J. This calculation relied on 
using spinor variables in the tree amplitudes composing the unitarity cuts. 

At loop level the application of four-dimensional helicity methods is much trickier than 
at tree level because of the use of dimensional regularization, which requires that the loop 
momenta be outside of four dimensions. This can drop rational pieces not captured by 
four-dimensional unitarity cuts. A purely four-dimensional method relies instead on on- 



shell recursion relations to reconstruct the rational pieces [10|, Ilia, 116] . Another approach 
is to abandon four-dimensional helicity in the cuts and use D- dimensional generalized cuts, 
which simultaneously give cut-constructible and rational terms 

geH. Rational terms 



arise through the interference of the divergences in loop integrals with the order e pieces of 
(4 — 2e)-dimensional loop momenta. A means for tracking these within the OPP framework 
in terms of a set of Feynman-like rules has been given in ref. [l9[. Work by Ellis, Giele, 
Kunszt and Melnikov combines D-dimensional generalized unitarity with the OPP method 



for numerical calculations |20H22|. A related approach due to Badger [23J associates the 
(D — 4)-dimensional pieces to masses to give the rational terms, and a numerical version of 
this has been developed and applied in refs. [H], 13 • 



Given the success of helicity methods when combined with generalized unitarity at push- 
ing back the frontiers [3] , it seems reasonable to try to find improved means of combining the 
two techniques. In this paper we bring D-dimensional unitarity and spinor helicity together: 
we work in D = 6 using a six-dimensional helicity formalism due to Cheung and O'Connell 
2^1 to calculate QCD amplitudes, including NLO corrections to Higgs processes, at one loop. 



We capture both the cut-constructible and rational pieces while enjoying the conveniences 
of a spinor-helicity formalism, and are in fact able to write the answer in terms of massless 
four- dimensional spinors through a decomposition of their six- dimensional counterparts. At 
the end of the calculation, we analytically continue the loop momenta to D = 4 — 2e dimen- 
sions and perform a state-sum reduction to reduce the spin states of six dimensions to match 
our desired regularization scheme. In this paper we focus on the four-dimensional helicity 
(FDH) scheme jB], [13, 26|| • Loop computations combining six- dimensional helicity and uni- 



tarity to obtain supersymmetric amplitudes have been carried out in refs. [27H29|. though 
for the cases treated there, there was no need to apply a state-sum reduction procedure as 
we do here. 

We illustrate the technique using various four-point QCD amplitudes as examples, in- 
cluding ones with external quarks. We also present a simple five-point example. As a more 
sophisticated example, we evaluate an amplitude involving the Higgs boson and three gluons. 
This process is mediated by a heavy-quark loop [30], so NLO involves a two-loop process. 
However, in the large top-quark mass limit, we can replace the top loop with an effective 



operator [3lH34|. thereby bringing our NLO computation down to one loop. In all cases we 
reproduce known results confirming the validity of the formalism. 

This paper is organized as follows: In section [Hj we review the six-dimensional spinor- 
helicity formalism of Cheung and O'Connell [25J. We cover the embedding of four- 
dimensional spinors into six-dimensional spinors, a crucial step to obtaining amplitudes 
in terms of four- dimensional objects. In section IIHI we outline our choice of integral ba- 
sis in 4 — 2e dimensions for the FDH scheme and identify coefficients that must be found 
to obtain the full loop amplitude. Then in section IIVI we review integral coefficient ex- 



traction for (4 — 2e)-dimensional generalized unitarity following Badger [23|. We give the 
appropriate loop momentum parametrization and coefficient formulas, leaving the details 
behind the procedure to an appendix. Section |V] brings us to sample calculations of the 
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one-loop four-gluon amplitude A\' (1~, 2 + , 3 + , 4 + ) for various internal-particle states. We 
find that a state-sum reduction is necessary to eliminate the extra-dimensional polariza- 
tion states of six dimensions compared to four dimensions. We also present our results for 
a maximally helicity-violating configuration for a gluon loop, as well as for the five-point 
amplitude A^\l + , 2 + , 3 + , 4 + , 5 + ). In section fVTI we discuss the case of external massless 
quarks, also subject to a similar state-sum reduction procedure. In section IVIII we apply 
our technique to a simple Higgs to partons NLO QCD correction process, a topic of great 
interest at the LHC. All Higgs plus four partons processes have been calculated analytically 
using four- dimensional cuts and various other methods for the rational pieces [35| . We show 
the utility of our method with the simple example A^'(H, l + ,2 + ,3 + ). In section IVlIII we 
give our concluding remarks. 



II. SIX-DIMENSIONAL HELICITY 

We begin with a summary of Cheung and O'Connell's six- dimensional helicity formalism 



25| , which gives us a spinor- helicity formalism to use in unitarity cuts where internal parti- 
cles must be kept in D ^ 4. After a brief review of four- dimensional spinor helicity, we make 
the analogous construction in six dimensions using solutions to the massless six- dimensional 
Dirac equation in a similar 'Weyl" basis. We then show how the six-dimensional spinors can 
be decomposed into four- dimensional spinors. In Section |V] we will explain how to correct 
for the state sums being in six dimensions instead of the dimensionally regularized values. 



A. Spinor helicity 

The four-dimensional spinor-helicity formalism has been widely used in scattering ampli- 
tude calculations, and a useful review can be found in ref. (36|. We use massless chiral and 
anti-chiral spinors in a Weyl basis to represent a light-like momentum p M as a bi-spinor, 

Paa = PtiKa = Kx^a- ( 2 -l) 

In addition to its momentum p^, a massless particle in four dimensions is labeled by its 
helicity. We define polarization vectors in terms of spinors X p = \p + ) and X p = \p~), which 
can be done as 

where is some arbitrary light-like reference momentum representing a gauge freedom. 
Following this formalism, we are able to write compact expressions for amplitudes in terms 
of Lorentz-invariant spinor products, 

e ap \ ip \ 3a = (ij) and e^Af \f = [ij], (2.3) 

where e 12 = €21 = 1. Here we are following the convention = Ipi ■ Pj = Sy. 

In six dimensions the Lorentz group is SO(5,l). We can then decompose the Dirac 
equation into chiral and anti-chiral pieces using a "Weyl" basis, 

P^AB*p a = 0> P^Kb, = 0, (2.4) 
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where {A, B, •••} are fundamental representation indices of the covering group, SU*(4). 
The sigma matrices oab and a AB are 4x4 antisymmetric matrices playing a part analogous 
to the Pauli matrices of four dimensions. Explicit forms and some useful relations can be 
found in Appendix A of ref. [25j. The chiral and anti-chiral Weyl spinors X Ba and Xsa each 
have two solutions labeled by the indices a = 1,2 and a = 1,2. These are indices of the 
little group SO(4), corresponding to SU(2)xSU(2). They can be raised and lowered using 
the matrices e ab and e ab , 

X a = e ab X b , A d = e<% (2.5) 

where once again e 12 = e 2 i = 1. 

Using these spinors we can write a six- dimensional light-like momentum vector in a bi- 
spinor representation similar to that of four dimensions, 

P AB = P^ AB = X Aa e ab X Bb = \p a )e ab (p% 

Pab = V^ab = XAae" b X Bb = \Pa]e db [p b \, (2.6) 

where we have adopted a bra-ket notation X Aa = \p a ), X A& = \p&]. We can also express the 
momentum vector directly as 

P" = ~\(p a W> b K b = -\[Pa\~^\p b Y\ (2.7) 

where there is a contraction of SU*(4) indices between the sigma matrix and the spinors. 
Lorentz- invariant spinor inner-products are also defined by contractions of SU*(4) indices, 

0. (2.8) 

Other important quantities showing up in amplitude calculations are the spinor contractions 
with the SU*(4)-invariant Levi-Civita tensor, 

(j a j b k c l d ) = e ABCD Xt a XfX^Xf d , 

ia.iiJ'.l.j. = e AB ° D 'XiAaX jBb X k Cc\ D d, ( 2 -9) 

and spinor strings, 

WAA ■ ■ -K^) = (^(P^MAM^ 3 ■ ■ ■ (P2n + l)A 2n+1 A 2n+2 (X j ) A ^ b , 

(i a \t x t 2 ■ --tjn] = {\) Ma {Pl)MA 2 {P2) MM ■ ■ ■ (P2n) A2 " A2n+1 (X 3 ) A2n+lb - (2-10) 

Finally we look at polarization vectors, which, as in four dimensions, can be written in 
terms of spinors. Following ref. (25| . 

= -^{p-KlftX^lP 4 ])" 1 = -^((p a l%])^[%I^N- (2.ii) 

Unlike the four-dimensional case, we cannot label these as simply + or — because gluons in 
six dimensions have four polarization states. The states are in fact labeled by SU(2)xSU(2) 
little-group indices. For the Weyl spinors X Aa and A^d then, the indices a and a label two 
helicity states respectively. We notice that the six- dimensional polarization vectors have a 
gauge freedom through the null reference momentum q, as is the case in four dimensions. 
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B. Decomposing six-dimensional spinors into four-dimensional ones 



We find it convenient to write six- dimensional spinors in terms of four-dimensional ones, 
allowing amplitudes to be expressed in terms of the more familiar four- dimensional spinors. 
From a four- dimensional perspective, we view six- dimensional null vectors as being massive. 
Making the associations, 

m = p 5 -ip 4l m = p 5 + ip 4 , (2.12) 

our six-dimensional massless condition is 

p 2 = p 2 — p\— p\ = p 2 — mm = 0, (2-13) 

where p denotes a momentum vector with only the first four components. This is the 
on-shell condition for a four- dimensional massive momentum. We can write a bi-spinor 
representation of massive momenta in terms of two pairs of four-dimensional spinors: A, A 

and fj,, p,, as 

Paa = KK + P PaP-w ( 2 - 14) 

where 

~_ /-/ _ m ~_ m , _ m ~/_ m 

p — KK — K K , K — — -, K — - -r, K — — -, K — - —, l / - i0 J 

(Xfi) [fi\\ (Xfj) [pX] 

which leads to a form for the six-dimensional spinors. Treating them as 4x2 matrices, we 
decompose them in terms of the above four-dimensional spinors as 27], |28|, |37 

The SU*(4) indices label the rows while the little group indices a and a take on the values 
1 and 2 to label the columns. The embedding is specific to the form of the a^ AB matrices as 
taken from ref. (2o| . 

If we take the momenta to be in the four-dimensional subspace with p 4,5 = 0, or equiva- 
lently m = m = 0, the spinors reduce to the four- dimensional forms, 

, / X a \ r ( X a 



A «=U" J' ^ M = \ —Xa J' ^ 

When we do unitarity calculations, we will keep the external particles in the four- dimensional 
subspace, so they will have this simpler form. The internal particles, on the other hand, will 
have the more complicated form of eq. (12.161) . 



C. Tree-level amplitude examples 

We now consider some tree-level amplitudes that we will need for unitarity cuts. Tree- 
level amplitudes in six dimensions have a remarkable chiral-conjugate structure. The color- 
ordered four-gluon tree amplitude is given by [25] 



^(1^,2^,3^,4^) = (U^lWj. (2.18) 
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This was found in ref. 25] by building u p f rom three-point amplitudes using the Britto- 
Cachazo-Feng-Witten recursion relations [38|. Notice that, unlike four-dimensional ampli- 
tudes, the helicity is not specified. The little group SU(2) xSU(2) connects all helicities 
together so this simple result holds for all helicity arrangements. For the two-chiral-quark 
two-gluon tree amplitude, we have [27] 



4 0) (1L 2^,32, 4*) = -— — (l a W d >[l,2 fe 3^], (2.19) 

-"12 323 

while for the two-real-scalar two-gluon tree amplitude, 

4 0) (1L,2^3 S ,4 S ) = --^(l a 2 6 3 e 3 e )[l a 2 b 3 e 3 e ]. (2.20) 

One can easily check that these reduce to the known four-dimensional results using the 
reduced spinors in eq. (12.171) . 



III. (4 - 2e)-DIMENSIONAL INTEGRAL BASIS 

A general one-loop amplitude can be written in terms of a basis of master integrals. We 
give our choice of basis here and leave the details to Appendix [Aj where relevant references 
can be found as well. Following ref. 23j , for n particles in four dimensions with all internal 
particles massless and in D = 4 — 2e, a one-loop amplitude can be written as 

4(1) _ ^ [ v r [0] T*~ 2e + V r [4] r 4 - 2e [/; 4 i + V r [0] T 4 ~ 2e 

K ' \ K 4 K 4 K 3 

+ E ^l^Vl + E ^1 J S + E ^Kjtrg\A) +0(e), (3.1) 

K 3 K 2 K 2 J 

where K r refers to the set of all ordered partitions of the external momenta into r distinct 
groups. The integrals are given by 



(2tt) 4 -^ l 2 (l - KtfQ -K x - K 2 ) 2 ■■■(l + K r , 



For simplicity in eq. ( 13. ip . we have written / 4_2e without the explicit /j, 2 dependence when 
f(/J> 2 ) = 1. The scalar integrals have been evaluated and can be found in refs. 0, 39]; 
the integrals with f(n 2 ) ^ 1 are straightforward to evaluate and are given in Appendix lAl 
Because all the integrals in eq. ( 13. ip are known, calculating one-loop amplitudes boils down 
to finding the coefficients in front of the integrals. 



IV. INTEGRAL COEFFICIENT EXTRACTION 

With up to box integrals in our basis, the largest number of cuts necessary for generalized 
unitarity in 4 — 2e dimensions is four. Quadruple cuts give us box integral coefficients, while 
three- and two-particle cuts give us, respectively, triangle and bubble integral coefficients. 
However, there is some entanglement depending on the cut. Two-particle cuts, for example, 
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are also contained in box and triangle integrals, and the bubble coefficient must be isolated 
through some extraction procedure. An efficient method for coefficient extraction in D = 4 
based on the analytic properties of the integrand is given by Forde js)]. Badger 23] generalizes 
this method to D = 4— 2e with an emphasis on rational terms obtained by treating the extra- 
dimensional components of momenta effectively as masses, in a form related to extraction 
of integral coefficients in the presence of masses [40] . We follow the construction of these 
references, providing the momentum solutions to the cut conditions and the final forms of 
the coefficients. Further details on the procedure can be found in Appendix 03] 



A. Box integral coefficients 



w 

h 

— , H — 

MS 



FIG. 1: A general quadruple cut. Loop momenta flow clockwise. 



For the quadruple cut (Fig. [TJ in 4 — 2e dimensions, the on-shell cut conditions are 

If = ll = ll = 11 = 0, (4.1) 

or equivalently, 

n=%=n=n=f, ^ 



where U are the momenta truncated to four dimensions and /i represents the (— 2e)- 
dimensional components. It is convenient to view fi 2 as a mass term and the internal 
momenta as four-dimensionally massive. (This is the view taken in ref. 23J and as we saw a 
useful one for six-dimensional helicity as well.) We parametrize the four- momentum l± using 
a construction first given in refs. 0,0, 4l[ and also applied in refs. [s], [23|, 40]. 



11 = -(a(Kl\^\Kl) + b(K'f\Y\K\-) + c(K\-\Y\K\~) + d(K\-\Y\K\-)), (4.3) 

where K\ 4 is the massless projection of one of the external legs in the direction of the other 
masslessly projected leg, 



1 7i4 - 



K 



b M _ 714(714-^4 ~ s * K i 
7i4 - 



S 



714 = K 1 -K i ± v^i • K,f - S { = Kf. (4.4) 

Solving the on-shell conditions, we find that l\ can be expanded as 

% = aKY + bKY + C -±(K\-\r\K\-) + 1u ^~^\ k\-\^\K\~) 

= fT-7r—(K b l -\^\Kl), (4.5) 
2c±7i4 



where 



Si{S4 + Ju) b _ S A (S 1 +714) _ -Ci ± a/c 2 - 4:C C 2 

7i4 - ' 7i4 - ' 2c 2 

c 2 = {K\-\t 2 \K{-), 

Cl = a{K\-\$ 2 \K\-) + b(K\-\$ 2 \K\-) -S 2 - 2K X ■ K 2 , 



(ab-£) (K\-\$ 2 \Kt). (4.6) 

In general there are two solutions to the on-shell conditions. It would appear initially that 
there are four, two each for 714 and c±, but it turns out that, 

h(lf 4 ,c + ) = fi(7r 4 ,c_), 

/i(7i + 4 ,c-) = l 1 (-Yu,c+). (4.7) 

In the case S\ = or £4 = 0, there is only one solution for 714 (but still two solutions to 
the on-shell conditions). To determine the full box coefficient, we must average over these 
solutions. 

Using these solutions to define spinors, we can calculate the four tree amplitudes associ- 
ated with the quadruple cut. As explained in Appendix [Bj the coefficients associated with 
our integral basis choice ( 13. II) are 

/ \ -i 4 , j.. 4 . / /" \ 

l/i 2 -K) ' 



cf = 1 X; A 1 A 2 A 3 A A (jr J 



= \ E [Inf^A^M/^U • (4.8) 

(7 

The sum is over the two solutions to the quadruple cut; the product A X A 2 A 3 A4 must be 
computed for each. To find Cf \ our cut-constructible piece, we set /i 2 = 0. The Inf term 
111 Cf, first given in ref. [42|], contains the information from the boundary of the /i contour 
integral, 

lim ([inf^A^A^ 2 ) - A x {ii)A 2 {ii)AMA A {ii)) = 0. (4.9) 

We expand around infinity and write it as a polynomial in /j, 2 , 

2 

[Inf M2 AiA 2 A 3 A 4 ](/i 2 ) = X ( 4 - 10 ) 
then restrict C4 to be the coefficient of the /z 4 term, explaining the notation of eq. (14. 8p . 
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B. Triangle integral coefficients 




To obtain the solution to the triple cut (Fig. [2]), we again parametrize our loop momentum 
in terms of adjacent projected external momenta, 



2 v ° " 1 1 ' 2i 7 



13 



where K\ and K\ are defined analogously to eq. (14. 4p . The coefficients a and b are defined 
in eq. (I4.6P and 713 is defined in eq. (I4.4j) (with K"g replacing K\ in both). We must average 
over the solutions for 7 13 , though in the case that S\ = or £3 = 0, there is only one 
solution. For a fixed value of 713, we must also average over the coefficients given by the 
conjugate solution, 

1*/ = aKf + bKY + -(K:\-\Y\Kl-) + llsab ~ ^ \ kI~\Y\K\-). (4.12) 



2 v 1 " 1 ' 2t~f 



13 



In both solutions the complex parameter t is free. 

Box integrals also contain triple cuts, so we must extract the triangle coefficients using 
the limiting behavior of the integrand. The coefficients therefore contain an Inf term that 
is a polynomial expansion in t, 



2n ^ ^ L t J- ^ ov l/JV /| ^^q^^q 

C f = [Inf, 2 [Inf^ 1 A 2 A 3 (/n](t)](/i 2 )| M2 ^ , (4.13) 

but only the order t° term is retained. The sum is over the solutions, including the conjugate- 
momentum solution, to the cut conditions. There may be either two or four solutions 
depending on the number of solutions n 7 for 713. In C3 ', fi 2 and t are both set to zero, while 
the expansion in Cg 2 ' is restricted to the coefficients of the fi 2 term (in addition to having 
t = 0). 
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FIG. 3: Double and triple cuts contributing to bubble coefficients. 



C. Bubble integral coefficients 

Two-particle cuts are contaminated by both boxes and triangles, so again an extraction 
procedure is necessary. Furthermore, triple cuts that share two of their cuts with the double 
cut contribute to tensor triangle integrals that reduce to scalar bubbles, so we must take 
them into account for the full bubble coefficient (Fig. [3]). The two-particle cut momentum 
solution has two free parameters t and y, 



7 A Zt^f 



where 



b/t 



7 



7 = 2(#i-x), 



(4.14) 



(4.15) 



and x is some arbitrary massless vector. We calculate the cut integrand A\A 2 by taking a 
product of the two on-shell tree amplitudes, which will make up part of our coefficient after 
an extraction procedure. 

For the triangle contribution to the bubble coefficient, we fix the parameter y to put 
another propagator on shell, 



y± 



B 1 ± ^Bf + AB B 2 



2B 2 



(4.16) 



where 



B 2 = S l ( X -\^s\K]r), 
B 1 = ^t(K b 1 -\tt 3 \K b 1 -) 



S 1 t{ X -\$^) + S l {x-\tz\K\-) 



B Q = *{t 2 (K\-\$z\x-) -f{x-\$z\K\-) +jtSs + tS x ( X -\tz\xr)- (4-17) 
We then calculate the triple-cut integrand for all triple cuts that share two cuts 
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with the original double cut. The bubble coefficients are then given by 



Ctri CTy 

cf = - i [Inf At2 [Inf 4 [Inf^ 1 A 2 ](y)](t)]( M 2 )| M2 ,^ ,^^K m 

" ^EE[ Inf M 2 [Inf^i^^3](t)](^)U 2 ,^T r (4.18) 

In the series expansion [Inf y Ai A2] (y) = Ylm=ofmy m , we make the replacements y m — >■ F m 
where 



\ >HK)> ^-K 1 - 2 ©- MHH' 

(4.19) 



These double-cut terms are also expanded in t and restricted to the t° term. In Cf \ fi 2 is 
set to zero, while in Ci, 2 ' there is an expansion in /z 2 , and is restricted to the order /i 2 
coefficient. The triple-cut terms have a sum over the two solutions for each triple cut and a 
sum over all possible triple cuts that share two cuts with the double cut. In the expansion 
in t, we replace P —> Tj where 

T = 0, 

S l { X -\t,\K\-) 



27A 

3^( X -|^ 3 |^') 2 
87W 

(\ 1jt 3 |^ -> 3 
48 7 3 A 3 



T 3 = ~ r^ A 3 ( 155 i 35 3 + 30ATi ■ tf 3 S?S 3 + ll(^i ■ K 3 ) 2 Sf + 4SfS 3 + lQ^SlA) , 

(4.20) 



and 



A = (K 1 - K 3 ) 2 - S 1 S 3 . (4.21) 



In a renormalizable gauge theory (such as QCD), double cuts can have terms up to order 
y 4 while triple cuts can have up to t 3 , so all necessary integrals are evaluated above. 1 Again 
see Appendix [B] for a more detailed explanation. 



V. COMPUTATION OF FOUR- AND FIVE-POINT GLUON AMPLITUDES 

With the formalism established, we turn to sample calculations of one-loop amplitudes in 
QCD, starting with ^^(l - , 2 + , 3 + , 4 + ) for a scalar loop, gluon loop and fermion loop. We 



1 As it turns out, y 3 and y terms come with factors of 1/i and 1/t 2 , respectively, so such terms drop out 
when the coefficient is restricted to the t° term, but we include the evaluations above so as not to break 
the procedural nature of the method. 
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use generalized unitarity to glue together tree amplitudes written in terms of six-dimensional 
spinors. In keeping with Cheung and O'Connell j25|, these six-dimensional spinors are 
written in terms of four-dimensional ones based on the loop-momentum solutions of section 
HVl The cuts then fall into a form that allows us to use our coefficient extraction techniques. 
Our methodology for the cuts follows that of ref. j27|. 

As mentioned earlier, internal particles are kept in six dimensions while external particles 
live in the four- dimensional subspace. This leads to an additional complication: increasing 
the number of spacetime dimensions also increases the number of spin eigenstates. We 
therefore do not expect the six-dimensional coefficients to exactly match the four- dimensional 
ones, and in fact we must perform a state-sum reduction to bring our result back to four 
dimensions. 

We find that the scalar-loop diagram does not require adjustment, nor does the fermion 
loop. For the gluon-loop diagram, however, we must subtract two factors of the scalar- 
loop diagram. We motivate these procedures by tracking how the dimension makes its way 
into the integral coefficients and demonstrate their effectiveness in the one-loop amplitude 
A { l\l-, 2+ 3+ 4+). We then include some additional results, but with fewer computational 
details. 



A. One-loop four-point cut conditions in D = 6 



2 3 2 3 2. .3 




(b) 



FIG. 4: Cuts for the four-point amplitude. We show triple and double cuts in the S23 
channel; in general we must also evaluate cuts in the s\2 channel. 



To compute the amplitudes using six- dimensional helicity, we use momentum solutions 
from section IIVI that fit with the coefficient extraction techniques. There we interpreted 
the (— 2e)-dimensional component \x as a mass and gave the four- dimensional components 
of the loop momenta the on-shell cut conditions, If = /i 2 . In section QI] we also expressed the 
extra-dimensional components of six- dimensional momenta in terms of mass parameters m 
and m given in eq. (12. 12ft . 

4 _ —m + fh 5 _m + m 

k ~ 2i , k ~ —3—, I 5 - 1 ) 

giving the on-shell cut conditions, If = mm. This enabled us to embed four- dimensional 
spinors in the six-dimensional ones. Making the association fi 2 «-» mm, we see that we have 



13 



a bi-spinor momentum representation that simultaneously allows us to use six- dimensional 
helicity and the coefficient extraction techniques. 

Using eq. (14. 5p . the momentum solution for the quadruple cut (Fig. |4aj) is 

(^)i = ~ (c±<4-|y|l-> - -4— (11714-)) , (5.2) 

where 



+ \ 12 ) L L . 4mmsi 3 



cf = -V-V 1±W1+ ■ (5.3) 

1 2(4 2) ^ V s 12 s 23 1 

The other loop momenta can be found through momentum conservation: l 2 = h — Pi, 
I3 — h — Pi — p 2 and I4 — li + p^, or by relabellings in eq. (15. 2p . It is easy to confirm that 
all are on shell. There are in total two solutions to the quadruple cut, and we must average 
over them at the end of our calculation. 

For the triple cut in Fig. [4bJ we use eq. (14. lip , 

n = U *<4-i7in - 7— aiTi^ , (5.4) 



2 \ t S14 

which solves the on-shell conditions l\ = l\ = l\ = mm. As before t parametrizes the 
remaining degree of freedom. Notice that requiring the remaining internal propagator be 
on shell (Z| = mm) recovers the box solution, i.e. t = cf. We will also need the conjugate 
solution for the cut in Fig. | 



IT = \ (t(l-|7l4-> - ^<4-|7in) • (5-5) 

Solutions for the other triple cuts can be constructed in the same manner. 

Finally for the double cut in Fig. HcJ we have two free parameters for l\ = l\ = mm. 
Choosing K\ = p\ + p^ and x = Pi m e Q- (14.141) . we have K\ = p$ and the momentum 
solution is 

y<4-[y|4-> + (1 - y)(l-|r|l-> + *<4-|y*|l-> + ~ y t J " (l-|7l4-) 

'(5.6) 

Setting y = 1 puts l\ on shell and, in fact, gives us the momentum solution to the triple cut 
in Fig. 



B. One-loop four-point solution for six-dimensional spinors 

With the momentum solutions above, it is straightforward to plug them into eq. (I2.16P 
and obtain six- dimensional spinors in terms of four- dimensional ones. We work explicitly 
with the double cut; spinors for the triple and quadruple cuts can be found by setting 
y — 1 and t = cf where appropriate. We start by expressing the solution in two-component 
notation, 

h = y A 4 A 4 + (1 - y) A1A1 + 1 A 4 A X + — ^— A1A4, 

T n \n \ \ . +\ \ , y(l-y) -mm/s u ~ 

h = -(1 - y) A4A4 - y A1A1 + t A4A1 H A1A4 . (5.7) 
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Making the associations, 

, ~ A 4 

fj, = X 1 , ^ = ~' 

\ — y y 

X h = —j—^i + A 4 , A r4 =tAi+yA 4 , A, 2 = --Ai + A 4 , Xj 2 = tAi - (1 - y)A 4 , 

Kl4= (41)' Kl4= [41]' " 14= (41)' Kl4= [41]' (5 ' 8) 
a comparison with eq. ( 12. 16ft gives us the spinors, 

i^h) a 
(AjjAd 
i^h) a 



-ki 4 Ai (1 -y)Xi/t + A 4 
tAi + yA 4 Ki 4 A 4 /t 

< 4 Ai _ (1 -y)Xi/t + X 4 
-tXi-yX 4 K'i 4 A 4 /i 

. i / -«i 4 Ai -yAi/t + A 4 
tAi-(l-y)A 4 K U X 4 /t 

,r x / «i 4 Ai -yXi/t + A 4 A 

(Ai2)Aa -l-tA 1 + (l-y)A 4 *' 14 A 4 /t J' (5 - 9) 

With these it is simple to work out the spinor products required for a cut. It should be 
noted that because we use an all outgoing convention, incoming momenta are labeled —p. 
For purely gluonic amplitudes, it is sufficient to add a factor of % to the negative momentum 
spinors to handle these cases. 2 As an example, the spinor products relevant to the gluon-loop 
contribution to ^(l", 2+, 3+, 4+) are 



(H2)i,2i,3i,(J4)i) = det 
= 0, 

((-/ 2 )i,2 1 ,3i,(/ 4 ) 2 ) = det 



— Ik,\aX\ — k 14 Ai 
itXi - z(l - y)X A A 2 A 3 tXi + yX 4 



((-/ 2 ) 2 ,2i,3 1 ,(/ 4 ) 1 ) = det 



((-/ 2 ) 2 ,2i,3 1 ,(/ 4 ) 2 ) = det 



-ikuXx (1 -y)Xi/t+ A 4 
itXi — i(l — y)X 4 A 2 A 3 Ki 4 A 4 /t 
m 14 (14)[2 3], 

-iyXi/t + iX^ — «i 4 Ai 

ikiiXi/t A 2 A 3 tXi+yXi 

-i« 14 (14)[2 3], 

-iyXt/t + iXi (1 - y)X x /t + A 4 

«i 4 A 4 /i A 2 A3 Ki 4 A 4 /t 



z~<14)[23], 



2 Cut fermions do not have this factor of i. We will generally define the spinor products with the factor of 
i in the cases of negative momenta and manually remove it when we cut fermions (see ref. fl~6| for more 
detail). 
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((-/ 4 )i, 4 1; 1 2 , (Z 2 )!> = 0, ((-khAi, h, (kh) = it(l 4) [4 1], 

((-/ 4 ) 2 , 4 1; 1 2 , (/ a )i) = -if (1 4) [4 1], ((-U) 2 , 4 1; 1 2 , (/ 2 ) 2 ) = 0. ^- iUJ 

For the bracket spinor products, we make the replacement k 14 — > k' u and change the sign 
where appropriate, 

[H 2 )i, 2 h 3 i; (Z 4 )i] = 0, [(-l 2 ) h 2 h 3 i; (h) 2 ] = -m' 14 (l 4) [2 3], 

[(-Z 2 ) 2 , 2 h 3i, (/ 4 )i] = zk' 14 (1 4) [2 3], [(-Z 2 ) 2 , 2i, 3 i; (Z 4 ) 2 ] = i±(l 4) [2 3], 



[(-/ 4 ) i5 4i, 1 2 , {l 2 )i\ = 0, [(-/ 4 )i, 4i, 1 2 , (Z 2 ) 2 ] = if (1 4) [4 1], 

[(-/ 4 ) 2 , 4 i; 1 2 , (/ 2 )i] = -if (1 4) [4 1], [(-/ 4 ) 2 , 4 b 1 2 , (/ 2 ) 2 ] = 0. 



(5.11) 



C. Scalar loop 

We are now ready to find the full amplitude for the scalar-loop contribution to 
^^(l - , 2 + , 3 + , 4 + ). As mentioned before, increasing the dimension of spacetime for the 
virtual particles in the loop increases the number of spin eigenstates, potentially causing a 
difference between coefficients in six-dimensions and those in the FDH scheme. However, 
scalar particles do not have spin, so we expect the coefficients here to match. This is indeed 
what we find. 

This fact can also be understood by examining the vertex factors in a one-loop Feynman 
diagram. The two-scalar gluon vertex is 

^(p-ffV. ( 5 - 12 ) 

where p and q are the momenta of the scalars. These momenta are contracted with the 
polarization vectors of the gluons, which, since they are external, are in four dimensions. 
Any extra-dimensional components are annihilated in these contractions, so whether we are 
in six or 4 — 2e dimensions, we get the same contribution from a given vertex. This argument 
can be extended to the two-scalar two-gluon vertex to handle all Feynman diagrams and 
show that there are no terms in the amplitude that are proportional to the dimension of the 
internal particles. 

Moving on to the computation, we begin with the quadruple cut (Fig. liaj) . The usual 
procedure is to sew four three-point amplitudes together, but working with three-point 
amplitudes in six dimensions can be complicated. It is, in fact, easier to multiply simpler 
four-point tree amplitudes by inverse propagators making them equivalent to the sum of 
products of two three-point amplitudes (Fig. [5]), given the cut conditions. The quadruple 
cut is then 

C 1234 = Yl M-ll, 1~ IDM-Il 2 + , IfiM-ll 3 + , ll)A 3 (-l s 4 , 4 + , ID 

states 

= H) 2 ( /2 - + PJ 2 M-l S 2, 2 + , 3 + , IfiM-li, 4+ , (5-13) 

states 

The two- real-scalar two-gluon tree amplitude was given in eq. (12.20 p . 

4 0) (1L2^,3 S ,4 S ) = -i-^-(l a 2 fe 3 e 3 e )[l a 2 b 3 e 3 e ], (5.14) 
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MA 



FIG. 5: The four-point quadruple cut. Two pairs of three-point amplitudes are grouped 
together (indicated by blue ovals) to form two easier-to-use four-point tree amplitudes. 
The cut propagators of the four-point tree amplitudes are canceled by multiplying by 
inverse propagators prior to imposing the cut conditions. 



so the quadruple cut becomes 

((-Z 2 ) a , 2 1; 3i, (-/ 2 ) a > [H 2 )d, 2 i; 3 is {-hf] (H 4 )&, 4i, 1 2 , {-k) b ) [(-Z 4 ) j, 4 ls 1 2 , (-h) b ] 



C 



1234 
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The relevant spinor products are 



(5.15) 



((- 


-h) 


i, 2i, 3i, i 




= <(■ 


-fa 


) 2 ,2i, 


3i, 


Ha 


) 2 ) = 


-k 14 (14)[23], 




((- 


-k) 


i,4i,l 2 ,( 




= (( 


-i* 


)2,4l 


la, 




) 2 ) = 


-t(14>[41], 




[(■ 


-k 


'i, 2i, 3j, 


(-k)'] 


= [(" 


-fe) 


2) 2j, 


3i, 




1 = 


4 4 (14)[2 3], 




[(■ 


-k 


)i,4 i5 l 2 , 


(-k) 1 ] 


= [(" 


-Z4) 


2A1, 


la, 




2 ] = 


-t(14)[41]. 


(5.16) 



Notice that we have not explicitly plugged in t = cf . We recommend that this be done 
only at the end of the calculation as the square-roots tend to fall out when we average over 
solutions. Multiplying the products out, we arrive at 



C 



1 



1234 



^ 14K ' 14 (14) 4 [4lH2 3]V. 



(5.17) 



'23 



Now we take our sign convention into account using eq. ( 14. 8 p and average over solutions 
that take Z 3 on shell, i.e. t = cf, arriving at 



C1234 

mm[2 4] 2 



('- 



•S12S23 / Si 2 S 23 



[12] (2 3) (3 4) [41] s 13 V 2s 



+ mm 



13 



(5.18) 



Associating mm — > fi 2 , we see that the cut is nicely in the form of a polynomial in fi 2 . We 
identify the coefficients, 



C 



c 



[4] 



[2 4] s 



Sl 2 S 23 



[12] (2 3) (3 4) [41] Sl3 



(5.19) 
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These are in agreement with the previously known result [5[. 

Next we examine the triple cut in Fig. |4b]to obtain triangle coefficients. Conveniently we 
have most of the cut for one of our momentum solutions from the box computation. We can 
use our quadruple cut in terms of the free parameter t and multiply by the l 3 propagator. 
This is the same as multiplying the three tree amplitudes in Fig. |4b] together, 

c 124 = A 3 (-i{, l- i$M-%, 2+ - 3+ > il)M-il, 4 + , ID 

= _ii_ %4/ ,; 4 (i 4 ) 4 [4 l] 2 [2 3)H 2 (t<4 2) [1 2] + s 12 - ^ (1 2) [4 2]) . (5.20) 

We also need the cut for the conjugate-momentum solution, 

IT = \ (f<l-|r|4-> - ^(4l7ll->) • (5.21) 

Instead of defining new spinors and spinor products, we recognize that the conjugate solution 
is obtained from the original solution through the substitution t — > —j^j- We can then make 
this same substitution in the cut, 

( - / 124 — ( - / 124|i— >— mm/tsi4 



-I 

S 



W; 4 (14) 4 [41] 2 [23] 2 ^^ft(12)[42] + Sl2 -^(42)[12]) \ (5.22) 

Finally we perform an expansion around t — > 00 as lined out in section IIVBI keeping only 
the order t° term, 

C 3 = -I {[Inf t (C m )]|^ + [Inf t (C* 24 )]|^ Q } 

1 f mm[2 4l 2 SioSoq 1 

J J r^ + 0L (5.23) 



2 1 [12] (2 3) (3 4) [41] S 2 3 
which gives us the triangle coefficients, 

r [o] _ n 

U 3;23 — U ' 



[2] _ i [2 4] 2 s 12 s 2 23 

° 3;23_ 2 [12] (2 3) (3 4) [41] sj 3 ' ' ' 



where "23" indicates that legs 2 and 3 have been "pinched" to form the triangle. 

We also need the triangle contributions from other triple-cut configurations. To pinch 
legs 3 and 4, we use the momentum solution, 

% = \ (Vl7l2-) - i^V| 7 11->) , (5.25) 

along with its conjugate solution. We define spinors and spinor products (many of which can 
be found through simple relabellings of our previous results) using this momentum solution 
and calculate cuts in a similar manner. The remaining triangle coefficients are then 

r [o] _ r [o] _ r [o] _ n 

°3;12 — °3;34 — °3;41 — u > 
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' 3 2 

°3 ; 12 - o m on /oov /„ 7\T< n I "T + 2 ~ S 23 

V s 13 S 13 

S 12 S 23 



,[2l 


i [2 4] 2 

L J 


3;12 ~~ 


2 [121(2 3) (3 4) [4 11 

L J \ / \ / L J 


-[2] 


z [2 4] 2 


3;34 _ 


2 [1 2] (2 3) (3 4) [41] 


-[2] 


i [2 4] 2 


3;41 — 


2 [12] (2 3) (3 4) [41] 



'13 

= ^ oTTo^buTTi ( f + 2^ - Sl2 ) • (5.26) 

'13 S 13 / 

Triangle integrals are dependent only on their kinematic invariant masses, so certain 
pairs of integrals corresponding to the above coefficients are in fact equal. For the one-mass 
triangle with invariant mass square S12 = (pi + P2) 2 = S34, 

= &f(v 2 )] = 4r 3 ?[f(n 2 )}. (5.27) 

It then makes sense to write only one coefficient C*3 ;Sl2 = C3.12 + C^a- In this basis our 
results are 



r [2] -i [24]2 523 (c -s ) 

° 3 ^- ~ [1 2] (2 3) (3 4) [4 1] S f 3 1513 12 J 



[2] [2 4]- s i2 

C ^ ~ '[12](23)(34)[41]4; (Sl3 " S23) - 



Finally we calculate the bubble integral coefficients using double and triple cuts. For this 
special case, dependence on y in the spinor products falls out and the products for these 
cuts are also given by eq. (15 . 16[) (this is of course not generally true). We can then take the 
quadruple cut (15.171) and multiply by propagators for Z x and Z 3 to obtain the two-particle 
cut, 

C24 =Ai(— 121 2 + , 3 + , l s ^)A±(—l\, 4 + , 1 , l s 2 



=Cl234 

1 



'2) 

i i 
II2 



til 



-—k 14 k' u (1 4) 4 [4 1} 2 [2 3] V(-(l - y)s u ) 

S 23 



ys 12 + (1 - y)sn + t (4 3) [3 1] + (^-A _ ^ (1 3) [3 4]) . (5.29) 

According to eq. (14.181) . the two-particle cut contribution to the bubble coefficient is 

-i [InU [Inf y A x A 2 \ {y)\ (t) \ t ^o, y ^Y m ■ (5.30) 
Expanding C24 around infinity in y, we find 

Inf y C 24 = 0, (5.31) 

so the contribution in this case is 0. 

We also need to examine three-particle cuts that share two on-shell conditions with our 
double cut, as explained in section HV CI The two possible configurations are shown in Fig. [61 
Fig. I6al makes no contribution due to our convenient choice of momentum parametrization. 
Any contribution must be proportional to 7\, T 2 and T 3 as given in eq. (14.201) . Recall 
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1 4 

(a) 



1 4 

(b) 



FIG. 6: Triple cuts contributing to the four-point amplitude bubble coefficient for the 

^-channel bubble integral. 



that for our momentum solution, K\ = p 4 , and in Fig. [6a], K 3 = p 4 . We then find that 
T± — T 2 — T3 = 0, and the contribution from this particular triple cut vanishes. 

The only bubble contribution then comes from Fig. l6b| where K3 = p%. To make use of 
eq. (I4.18p . we need to find the solutions to the parameter y that take the third leg 1% on 
shell. These are given by eq. ( 14.161) . 



where 



B x ± JBf + 450^2 
V± = V 2 ^ — , (5-32) 



B 2 = S23 <13)[34], 

B 1 = S23*(S12-S1 3 )+S23(13)[34], 

B = s 23 t 2 (4 3)[3 1] - fi 2 (l 3) [3 4] + s l3 s 23 t. (5.33) 



Our cut is given by 



c 2M = A 4 (-i\, 4 + , r, i' 2 )M-i a 2, 2+ , i s 3 )M-i s s, 3 + , il) 

= C123472 

= -^ Kl4K ' 14 (l4) 4 [41] 2 [2 3] 2 t 2 (-(l -y)s u )-\ (5.34) 
Plugging in the solutions for y and expanding in t, we find 

iEN«0 = l "'"<«34W-s 12 t + 0( ^ (M5) 

Lastly we need to take t — > T\, where T% is given by eq. f)4.20p . 

r , = -2<14M. (5.36) 

s 23 
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We then have the full coefficient expression, 

C 2 = -tilnUlIniyA^ymit^^ - \ £ J2[^tA 1 A 2 A 3 ](t)\p_ Tj 



2 

Ctri °~y 



= i [12](23>(34)[41]i^ (Sl3 - Sl2) ' (5 - 37) 
resulting in the coefficients, 

°2;23 — U ' 

°^ = ^12](2 3K3 4)[41] Sl 3 1 S 2 23 (Sl3 ~ Sl2) - (5 - 38) 

Through a similar computation, we find the coefficients for an invariant mass in the Si 2 
channel, 

Ly 2;12 — U ) 

2 ' 12 [12](2 3)(3 4)[41] Sl2Sl3 l 23; 1 ; 



All coefficients are in agreement with the previously known result [5 



D. Gluon loop 

Next we look at 2 + , 3 + , 4 + ) with a gluon loop. Unlike the scalar-loop amplitude, 

the gluon-loop amplitude requires a state-sum reduction to reduce from six dimensions to the 
FDH scheme. Gluons in six dimensions have two additional polarization states compared to 
four dimensions. The extra contribution of these states manifests itself in a term proportional 
to the dimension of spacetime, which is then the term that we must reduce. We will find 
that we must subtract twice the contribution of the scalar-loop diagram to reduce from six 



dimensions to the FDH scheme. The reasoning behind this can be found in ref. |2lj. We 
summarize the argument before presenting the results of A^\l~ , 2 + , 3 + , 4 + ). 

The dependence of a one-loop coefficient on the dimensionality of spacetime is at most 
linear. Such terms arise from a closed loop of contracted metric tensors and/or gamma 
matrices from vertices and propagators making up the loop. We can therefore separate 
coefficients into pieces dependent on D and pieces independent of it, 

= d + DC £>. (5.40) 



As the authors of ref. [21| point out, this can be used to find the coefficients in the FDH 
scheme by finding the coefficients in two dimensions D > 4: Di and D 2 (doing the full 
calculation in D = 4 of course loses the non-cut-constructible terms), then solving for Cj 
and Cd, 

<7(£>2) _ /7CD1) 
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The coefficient in the FDH scheme is then 
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C tuti = ^ r±r — + 4 r _ — . (5.42) 

D 2 -D 1 D 2 -D 1 V ; 



With this fact in mind, the authors of ref. [21| study the gluon-loop amplitudes in five 
and six dimensions and show that they differ in a rather simple way. Suppose we calculate 
in D = 5 with five- dimensional loop momentum. We will have three polarization states of 
the internal gluons, and we will obtain coefficients for five dimensions. Now suppose we 
upgrade to D = 6 but leave our loop momentum in five dimensions. We view this as being 
the same loop momentum as in the D = 5 case; the only difference is now we have added an 
additional polarization state of the gluon, which can be thought of as being along the sixth 
dimension of six- dimensional space. 

It should be noted that the restriction on the loop momentum does not prevent us from 
obtaining the correct D = 6 coefficient; problems arise only when we take it down to 
four dimensions. This can be understood through a fact mentioned in Appendix |Aj the 
loop momentum components beyond four dimensions only enter coefficients in a particular 
combination through contractions of the loop momentum with itself, 

p = I 2 - p, (5.43) 

where I and / are the four- and (D— 4)-dimensional pieces of the loop momentum respectively. 
As long as we do not set this to 0, we find the full coefficient (for a more in depth discussion 
see ref. [HI). 

So, the contribution to the D = 6 coefficient due to the extra polarization state (with 
polarization along the sixth dimension) is exactly the difference between the D = 5 and 
D = 6 coefficients. To see what this difference is, we go to the Feynman rules. 

A relevant three-gluon vertex, with polarizations taken into account, is 

A=(rT(k - l 2 ) p + V up (k - hT + rT{h ~ k) v )e kli e hv e hp , (5.44) 

where k is an outgoing external momentum and l\ and l 2 are outgoing internal loop momenta. 
We wish to examine the contribution from the polarization in the sixth dimension, so we 
look at only those components of and ei 2 . When these are contracted with e&, k, l\ or l 2 , 
they vanish since none of those vectors contain a sixth-dimensional component. The only 
non- vanishing contraction is that between the sixth-dimensional components of t\ x and q 2 . 
Our vertex then becomes 

e/i,6 ■ Q 2 ,6 efc ■ (h - h) = • (h - h), (5.45) 

up to a sign depending on the normalization of the polarization vectors. This is the two- 
real-scalar gluon vertex with polarization of the gluon included. Therefore the contribution 
of the polarization in the sixth dimension is equal to that of a scalar particle, and the D = 5 
and D = 6 coefficients differ by a coefficient found using cuts where the internal gluon loop 
has been replaced by a scalar loop, 
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For our FDH coefficient then, we find 



C FDH = C {6) - 2C scalar . (5.47) 

Coefficients in the FDH scheme can be found from six-dimensional coefficients by subtracting 
twice the contribution of a scalar-loop diagram. (Here we only studied the three-gluon 
vertex. Similar arguments can be used for a four-gluon vertex; in fact, it is true that any 
n-gluon tree amplitude with two gluons polarized alon g th e sixth dimension is equivalent 



to a tree amplitude with n — 2 gluons and two scalars [21].) This leads us to the physical 
interpretation that the polarization states in six dimensions consist of the polarization states 
of four dimensions plus two scalars, but that they may be mixed in some nontrivial way. We 
will not necessarily be able to identify the state a = 1, a = 1 in six dimensions, for example, 
with a state in four dimensions. 

We now turn to the coefficients of A^\l~, 2 + , 3 + , 4 + ) with a gluon loop and show that 
our state-sum reduction procedure works. Using the four-gluon tree amplitude (I2.18p . the 
quadruple cut is 



c-1234 = J2 H) 2 (^ - ^) 2 (/ 2 + pi) 2 a,(-ii 2 + , 3 + , ifiM-il 4+ > 1- > '!) 

states 

((-h)a, 2i, 3i, (U) b ) [(-l 2 )a, 2i, 3i, (h) b }((-h) b , 4 1; 1 2 , {l 2 ) a ) \{-U)\ 4 i; 1 5 , [l 2 f] 



6 23 



(5.48) 



The spinor products required for this calculation are in eqs. (I5.10j) and (15.111) . In a manner 
similar to the scalar-loop calculation, we find all of the integral coefficients for the gluon-loop 
amplitude. All scalar integral coefficients vanish, and we have 

c^' 6D = u-. ^ 24]2 512523 



[12] (2 3) (3 4) [41] Sl3 
r [2],6S [2 4] 2 s 2 23 

- 4 '[12](23)(34)[41]4; (Sl3 " Sl2) 

" 4! pppp5 l 23) 



n [2],6D_ A . [2 4] 2 S 2 3 / \ 

°^ - 4Z [12](23)(34)[41] Sl2Sl3 (Sl3 " S23) ' 

= 4z r w [ \f xr , 512 (s l3 - s 12 ). (5.49) 
2 ' 23 [12](2 3)(3 4)[41] Sl3S23 l ^ 1 ; 

Rational terms for a gluon-loop amplitude are equal to twice the rational terms for a scalar- 
loop amplitude [23] , so it is not surprising that the six- dimensional coefficients for the gluon 
loop are 4 times those for the scalar loop (this is not true for cut-constructible pieces). 
Subtracting twice the scalar- loop coefficients in section IV CI from eq. (I5.49P according to our 
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prescription, we arrive at our coefficients in the FDH scheme, 



Cf = 2* 



[2 4] 2 s 12 s 23 
[12] (2 3) (3 4) [41] Sl3 



°3; S i2 - r-i m/ nn > /o TTH n ~ l g 13 ~ Sl2j 



[12](2 3)(3 4)[41] S ? 3 

d 21 -2i [24]2 a "L- S 1 

( " 2 '[12](2 3)(3 4)[41]4; (Sl3 " S23) ' 



°^ ~ 2 ^12](2 3K3 4)[41] Sl S 2 2 S 3 13 (Sl3 " S23) ' 

°^ = 2 \l2}(2%%A)[Al}sZj SU ~ Sl2) ' (5 - 50) 

E. Fermion loop 

The final particle that we can have in the loop is a fermion. Fermions in six dimensions 
also have additional spin states compared to their four- dimensional counterparts, but we 
compute using six- dimensional chiral quarks with only two states, eliminating the need for 
a state-sum reduction. If non-chiral quarks were used, the loop amplitude in six dimensions 
would differ from that in four dimensions by a factor of two. This can also be understood 
through the Feynman rules. A quark loop gives rise to a trace over gamma matrices. Gamma 
matrices in four dimensions are 4x4 while in six dimensions they are 8x8. Where in four 
dimensions we have a factor of 4 then, we will have a factor of 8 in six dimensions. However, 
when we calculate using six- dimensional chiral quarks, gamma matrices are replaced by 
Cheung and O'Connell's sigma matrices, which are 4x4 and therefore contribute the same 
factor as in four dimensions. 

The two-chiral-quark two-gluon tree amplitude was given in eq. f l2.19p . Using it for the 
quadruple cut, we have 

C1234 = H) 2 (fe - ^) 2 (/ 2 + Pi) 2 A,(-ll 2 + , 3+ l\)A A {-ll 4 + , 1" l q 2 ) 

states 

_ (Hgjo, 2 U 3 U (k) b ) [Hgja, 2 h 3 i; {-l 2 f} ((-U)\ 4 1; 1 2 , (Z 2 )°) [(-k) d , 4 i; 1 2 , 

(5.51) 
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Following our earlier calculation methods 3 , we arrive at 

[2 4] 2 s 12 s 23 



[12] (2 3) (3 4) [41] a 



13 



- 2i [24]2 s ) 

°3; S i2 - ~ Zl T\ ol/oqWq AW A ll 72~l Sl3 ~ Sl2j 



[12](2 3)(3 4)[41] S 



13 



- 2i 12 4]2 s ) 

C ^ - ^[12](23)(34)[41]^ (Sl3 " S23) ' 



" - 2 ^[12](2 3)(3 4)[41] Sl2Sl3 (Sl3 " S23) ' 
°^ = " 2 ^12](2 3K3 4)[41] Sl S 3 1 S 2 23 (Sl3 " Sl2) ' (5 ' 52) 



where all other coefficients vanish. 



F. ^ 1) (l-,2 _ ,3 + ,4 + ) with gluon loop 

Because there are no cut-constructible pieces, A^\l~, 2 + , 3 + , 4 + ) with a gluon loop is 
just twice that of the amplitude with a scalar loop. Here we demonstrate the effectiveness of 
six-dimensional helicity in computing an amplitude with both cut-constructible and rational 
pieces. Specifically we look at A^\l~, 2", 3+ 4 + ). 

The six-dimensional coefficients for the gluon-loop amplitude are 

r [0],6D _ WO) r [4],6D _ dA (0)s 23 

4 — — /i 4 S i2 S 23 , 4 — — 4A 4 — , 



where 



r [0],6D _ n W0],6L> _ n 

°3;si 2 - U > °3; S2 3 ~ U ' 

G 3;,i 2 - U > G 3; S23 ~ U ' (5-53) 

^[01,60 _ n r [0],6D _ _io WO) 

2;12 — ' 2;23 — 3 4 ' 

^[21,60 _ n r [2],HD _ A W0) ( 2 1 , l" 

2;12 -U, <-2;23 ~ 4/i 4 I _ 3^ + HI 



(0) _ (12) 4 
^ 4 "'(12) (2 3) (3 4) (41)- (5 ' 54j 



For the scalar-loop amplitude, we have 



r»M _ _4(o)s23 r»[o] _ 1 4(0) 

O4 — /1 4 si2 , <-^2;23 ~~ 6 4 ? 



Wo. 



[2] 



2;23 



4 \ 3 S23 S12 J ' 



(5.55) 



where all other coefficients vanish. Subtracting twice the scalar-loop coefficients from the 
gluon-loop coefficients, we arrive at our coefficients in the FDH scheme, 

G 4 - -A A s 12 s 23 , G 4 - -2A 4 — , 

r [0] _ 114(0) r [2] _ 94 (0)/ 2 1 , i\ 5b J 

°2;23 - —^4 > ^2;23 ~ Z/1 4 ^ 3 i^ + ) > 



3 As mentioned in the previous footnote, we must manually remove a factor of i for each cut fermion. With 
two cut fermions here, we manually insert a negative sign in the computation. 
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in agreement with the previously known result [5| 



G. aJ 1) (1+,2+,3 + ,4 + ,5+) with gluon loop 



1 5 
FIG. 7: A possible quadruple cut for a five-point amplitude. 

As a five-point example we look at A 4 \l + , 2 + , 3 + , 4 + , 5 + ) with a gluon loop. For one- 
loop four-point amplitudes, we multiplied two four-point tree amplitudes together along 
with inverse propagators to obtain the quadruple cut. In a similar fashion, we multiply a 
four-point tree and a five-point tree together along with inverse propagators to obtain the 
quadruple cut here. A simple form for the five-point tree amplitude is given in section VII 
of ref. |25|. It was found through Britto-Cachazo-Feng-Witten recursion and repeated use 
of the five-point Schouten identity. However, we find it convenient to use an intermediate 
form before the use of the Schouten identity, 

^5 { L ad^ bb ^cc^ dd ^°ee) 



lT4r3f 4 



5,]) 



S12S23S34S45S15 

x |(l a 2 6 3 c 4 d )[l (i 2j ) 3e4 (i ](si4(5 e |^ 3 ^ 4 |5e] - si 3 (5 e |^ 4 ^ 4 |5j + si 2 (5 e 

- (l a 2 6 3 c 4 d )[l i 2j3 d 5,]( Sl 4<5 e |^ 1 ^ B |4 <i ] - s 13 (5 e |^ 4 ^ 5 |4 d ] + s^p^J^A,}) 
+ (l a 2 6 3 c 5 e )[l d 2^4 (i ](s 14 (4 d |^ l ^ 3 |5,] - s 13 (4 d |^^^j5,] + suiMpJJ^i]) 
+ s 15 (l 2 6 3 c 5 e )[l d 2j 1 3a5e](4 d |^ 1 ^ 3 |4 (i ] - s 45 (2 b 3Ad5e}[2 h 3^e}{la\pJJ 3 f 2 \U 

+ s 45 (2 6 3 c 4 d 5 e )[U2j ) 3a4 (| ](l ct |^ 3 ^ 1 |5e] + s 45 (2 fc 3 c 4 d 5 e )[l (i 2 i 3 £ 5e](la|^ 1 ^3|4j 

- s 45 (l a 2 b 3Ad)[2 b 3dM(5e\pJJ 3 p 2 \la} -s 45 (la2 fe 3 c 5 e )[2fe3,4 d 5e](4 d |^ 2 ^ 5 |l (i ]} . 

(5.57) 

The five-point tree amplitude used to find the quadruple cut will have three adjacent external 
gluons of positive helicity. Choosing these to be particles 2, 3 and 4, the form in eq. ( I5.57P 
leaves us with only one non-vanishing term, 



a(°)(-]9 o9,4D ,g,AD r 

^5 \ L aai z ii ' °ii '^li ' d 



ee) 



-(l a 2 1 3i5 e )[l d 2 i 3 i 5e](4 1 |^ 1 ^ 3 |4 i ], (5.58) 



S12S23S34S45 

where we have indicated that the external particles reside in the four- dimensional subspace. 
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There are five possible quadruple cuts for a five-point amplitude. We will calculate the 
cut in Fig. [7] explicitly and obtain the rest through cyclic symmetry. The quadruple cut 
momentum solution is 



T± _ r ±\ \ 1 mm r 
cf S15 



(5.59) 



where 



(12) 
2(5 2) 



1± W1 + 



4mms25 \ 



S12S23 



(5.60) 



The cut is given by 

CW = X)(-0 2 ^5Hf,2 + ,3 + ,4 + ,i2)^4(-lf,5 + ,l + ,Zf) 

states 

1 



S23 s 34Sl5(^4 +P4) 2 



((-Z 2 ) a , 2 1; 3 1; (/ 4 )6)[(-/ 2 )d, 2 i; 3 4 , (h) b ] (4 1 |/ 4 (-/ 2 )^ 3 |4i] 



x ((-Z 4 ) b , 5 1? l l5 (Z 2 ) a > [(-Z 4 ) 6 , 5i, li, (Z 2 ) d ]. 



(5.61) 



This does not fall conveniently into a polynomial in /1 2 , so we must use the extraction 
formulas in eq. (14. 8p . which give 



|,6D 
4;34 

[4],6D 



0. 



L/ 4;34 



4i 



tr + [3452] 



(12)(23)(34)(45)(51) tr 5 [1234] 



S 12^15, 



(5.62) 



where tr 5 [A] = tr^A] and tr+L4] = §tr[(l + 75) A]. After subtracting twice the scalar 
contribution, we arrive at 



W;34 — U ! 



c 1 ;. 



2i 



tr+[3452] 



4;34 



;i 2) (2 3) (3 4) (4 5) (5 1) tr 5 [1234]" Sl2Sl5 ' (5 ' 63) 

in agreement with the previously known results given in refs. (El, HH, 43]. All triangle and 
bubble coefficients vanish, and the remaining box coefficients can be found through cyclic 
permutations. 



VI. EXTERNAL FERMION AMPLITUDES 

Next we turn to the case of external fermions. The leading-color contribution to a one- 



loop amplitude with two external fermions can be decomposed as [44 



K = A L n - j^A* + ^A^M + ^, [0 ] (eu) 



c 



Only amplitudes with adjacent external fermions contribute at leading color. A^ is the 
"left-turning" fermion piece where the fermion follows the shortest path through the loop 
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FIG. 8: External fermion box diagrams for (a) the right-turning fermion and (b) the 

left-turning fermion. 



i 
i 

/"""X 

(b) 



FIG. 9: External fermion box diagrams with internal gluons replaced by scalars. 

(see Fig. l8b|) . while for the fermion follows the longest path (Fig. [Ha]). A^ 1 ^ and An'^ 
have, respectively, a fermion and scalar loop and are therefore proportional to the number 
of each. We focus on the left- and right-turning pieces and calculate the cuts in Fig. [HJ (For 
our sample helicity arrangement to follow, An'' 1 ^ and An'^ happen to vanish.) 

Once again the extra spin states of six dimensions must be treated with a state-sum 
reduction procedure. Fortunately our findings in section [V] give us good physical intuition 
hinting at the appropriate procedure. Fermions have additional states, but since we calculate 
with a six-dimensional chiral quark with only two states, no reduction procedure is necessary. 
For gluons, we interpret the extra-dimensional polarization states as scalars and subtract 
twice the contribution of an amplitude with a scalar loop replacing the gluon loop. The 
natural generalization then is to replace any internal gluons by scalars, even if they do not 
comprise the entire loop, and subtract twice the contribution of this diagram from a diagram 
with internal gluons. For example, cuts from the diagrams in Fig. [9] should be subtracted 
from their corresponding diagrams in Fig. [8] to reduce to the FDH scheme. 

For a more rigorous motivation, we again appeal to the Feynman rules. The dimen- 
sional dependence for Feynman loop diagrams involving external fermions arises through a 
contraction of gamma matrices. If we isolate the terms with dimensional dependence, we 
find them to differ only by a factor of D from equivalent diagrams with all internal gluons 
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replaced by scalars. We can therefore summarize our procedure as 



C FDH = C (6) - 2C scalar . (6.2) 
A. Right-turning A%\l+, 2", 3^ , 4+) 

As an example we calculate A4 (1 + , 2~, 3^, 4+) for both the right-turning and left-turning 
configurations. For the right-turning, we require the two- chiral- quark two-gluon tree am- 
plitude ( I2.19P to sew four-point amplitudes across the s 2 3 channel. We also need the four- 
chiral-quark tree amplitude with a gluon propagator 4 , 

4° );9 (M,3*,4*) = -z^Ml a 2 h 3 c 4 d >, (6.3) 

for sewing across the channel. Using these in our normal procedure, we find the six- 
dimensional coefficients, 

r [0),6D _ i 4 (0)£?2£23 W4],6D _ n 

^[01,60 = 1^(0) / _ s 12 s 23 \ C 12],6D _ A (0) f 2 sis _ £23 

3;12 2 4 ^ 12 s 13 J > 3;12 4 ^ "T" S23 Sl3 

n [Q]fiD _ 1.(0) n [2],6D _ AO) 

°3;23 — 2^4 ai , ) °3;23 ~ ^4 s , 3 > 

,|., ;; \ rt[2],6D _ i(0)/, , 

S13 l S23 .4 

°3;34 - ~2 A 4 (^12 + J , <^3;34 ~ A 4 ^ + + ^ > V ; 

^[0],6D _ i.(0) 

Sl2^23 

r,[2],6D _ AO) 
°3;41 - 2^4 Sl3 ' °3;41 ~ ^4 S13 > 



,6D _ 1 „ (0) , la , M W2],6D _ _A0) , 

«13 ' °3;41 - ^4 - 813i 

,loj,ei» _ _ 4 (0) ^[2]^ 

2;12 — ^4 ? Ly 2;12 

,[0],6D _ n r [2]fij 

2;23 ~~ U ' °2;23 



where 



- i (2 3)3(2 4) (6 5) 

4 ~ (12) (2 3) (3 4) (41)" 

To perform the state-sum reduction, our required tree amplitudes are the two-chiral-quark 
two-anti-chiral-quark amplitude with a scalar propagator (with adjacent chiral quarks) and 
the chiral-quark anti-chiral-quark scalar gluon amplitude, 

4 0);s (l«,2^,3l,4p = -^-(l a |4 d ](2 fc |3 c ], 

4 0) (M,3 s ,4p = -— J (l a 2 6 3 c 3 c )[l,3 c 3 c 4 d ], (6.6) 

4V2Si 2 S23 



4 If we wanted to compute a four-dimensional tree amplitude using this, we would need to reduce the state- 
sum by subtracting two times the diagrams with a scalar propagator replacing the gluon propagator. The 
six-dimensional chirality of the external quarks is not important in four dimensions since each chirality 
separately carries both helicity states of the four-dimensional quark. 
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as computed from Feynman diagrams. The coefficients for diagrams with internal gluons 
replaced by scalars are 



fnl 

cf = 


= o, 


cf = 


= o, 


°3;12 


= 0, 


°3;12 


= |4 0) ( 


r [o] 


= 0, 


r m 

°3;23 


— 4^4 


r [o] 

°3;34 


= 0, 


r m 

°3;34 


— 4^4 


r [0] 
°3;41 


= 0, 


r [2] 
°3;41 


— 4^4 


r»[°l 

°2;12 




r m 

°2;12 


= o, 


Ly 2;23 


= 0, 


r m 

Ly 2;23 


= 0. 



£12 

i aa 4. S23^\ (6.7) 

^ S23 ^ S13 / ' V ' 



£12 
S13 ' 



After utilizing our state-sum reduction procedure and combining terms with invariants in 
the same channel, we finally have, 

r [0] _ 1 4(0) gfgjgg ^[4] _ n 

W — 2^4 si3 ' °4 — U > 

- a(0) s 12 S23 -_/l(0)£23 

°3;si2 — ^4 S13 ) °3;si2 — ^4 S13 5 
^[0] _ „(0) S12S23 

°3;s 2 3 - A 4 ^3;s 23 ~ ~ A 4 J^) 

r»[°] _ _3 4 (») rrp] _ n 

L/ 2;12 ~~ 2 4 ) °2;12 ~~ U 5 

WO] _ n W2] _ n 

2;23 ~~ u ' °2;23 ~~ U ' 



in agreement with the previously known result [45] ■ 



B. Left-turning 2", 3^ ,4+) 

For the left-turning diagram (Fig. l8b]) . the only new tree amplitude that we require is 
the two-chiral-quark two-scalar amplitude, 



A 



(o), 



is 9S oq Al 



H g 13 - g23j 
4Si 2 S 2 3 



;iai a 3 c 4 d ). 



(6.9) 



Proceeding as usual, we find for our coefficients, 



C\ 

a 



— /1 4 (Si 2 S 2 3H 



3;si2 
,[0] 
3;s23 



A 



(0). 
4 

-A 



2 si 



■12^23 
«13 



(0) 



2;12 
2;23 



Sl2^23 
«13 : 



2 4 J 



0. 



c. 



[2] = 
3;si2 
W2] 

°3;s 2 3 - 
W2] 
°2;12 

W 2 1 
°2;23 



/l(O) S 23 

4(0)«ia 



0. 



(6.10) 



once again in agreement with ref. [45 



VII. HIGGS TO PARTONS PROCESSES: a£\h, 1 + , 2+, 3+) 

An important process for the Higgs boson search is production via the strong interaction. 



In the Standard Model, the Higgs couples to two gluons via a heavy-quark loop [30|. The 
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dominant contribution is from the top quark; contributions from other quarks are suppressed 
by at least a factor of 0(ml/mf) where m& is the mass of the bottom quark and rat is the 
mass of the top. 

Computations that keep explicit dependence on m t can be difficult, especially at NLO 
since they involve two-loop amplitudes. Fortunately in the large-m 4 limit, where the mass 
of the Higgs m# is smaller than the threshold for top-quark pair production, m# < 2m t , 



computations are greatly simplified by letting mt — > oo [31H33J. Accurate approximations 
(to within 10% of the Higgs-mass range at the LHC) of NLO QCD corrections can be found 
by finding the NLO QCD correction factor, K NLO = a NLO /a LO in the limit m t — >■ oo, then 
multiplying by the exact leading-order cross section |46|. The approximation is valid for 
multi-jet processes provided the transverse momentum of each jet is less than mt (47| . 

Taking the limit m t — > oo allows us to use an effective Lagrangian to express the coupling 
of the Higgs to gluons (33| , 

C^ = ^HtiG, u G^. (7.1) 

This reduces two-loop calculations at NLO to one-loop calculations. The Higgs two-gluon 
color-ordered vertex is 

-ZiW^kx-kt-h!?^), (7.2) 
which gives the six-dimensional helicity tree amplitude, 

4 0) (^,l«A,2 64 )=i(l«|2 i ]<2 6 |l a ]. (7.3) 
The Higgs three-gluon vertex factor is 

-^(rj^ih - k 2 y z +r]^ 3 (k 2 - h)^ 1 +rj» 3fll (k 3 - k^f 2 ), (7.4) 
which, when combined with other relevant Feynman diagrams, gives the tree amplitude, 



Af{E, lad, 2^) = 2V{(si 2 m| + Sl 3S 2 3)(l a |4|2 fe )[ld|^ 3 |2 ](3 c |^i|3c] 

Si2S 23 S 13 

- s 2 3(si 3 m 2 H + s 12 s 2 3)(la|3c](3 c |l a ](2 6 |^ 1 4|2b] 
+ si 3 (s 23 m^ + si 2 si 3 )(2 6 |3c](3 c |2j ) ](l Q |^2^3l 1 d] 
- Sl2 (l a |4fe|l a ](2 6 |^ 3 |2j ) ](3 c |^ 1 |3 6 ]}. (7.5) 



We use these to calculate A 3 >(H, 1+, 2+, 3+), first given in ref. [481 ]. 



A. Momentum solutions 

Momentum solutions can be a bit more complicated with a massive leg. We give explicit 
forms in line with those in section IIVI For the box in Fig. [TUl we use 

— ~ 1 TflTfl ~ 

k = c±AiA 2 A 2 Ai, (7.6) 

C± S12 
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FIG. 10: A possible box diagram for A 3 (H, 1 + , 2 + , 3 + ) using the effective coupling of the 
Higgs to two gluons. We denote this cut with ordering "123H". 



where 

The triple-cut solution for cutting propagators li, l 2 and l 3 can be found by letting c± — > t. 
To cut li, l 3 and / 4 , we use 

U = — — A3A3 + tXsX^H — A^A 3 , (7.8) 

13H tl3H 



where 



13H 



This allows us to set 



7 3// = 2k H ■ h 

= Sl3 - s 2 3- (7.9) 

A TfX^ 

[i = A#, [l = — , A = A3, A = — —A3 + tX H , 

t 13H 

m _ rh f m , m 

K3H= £xYy Kaa= w K3H = W^Y " 3H = ^xYY ( 0) 

Making similar associations for li and Z3 through the use of momentum conservation allows 
us to find the necessary spinor products. We find the cut for the conjugate momentum 
solution by taking t -> --™zl. 
To cut l\ and Z 3 , we use 

l x = yXtXx + (1 - y)X 2 X 2 + tAiA 2 + y{l - y) ~ mrh/Su X 2 ~X U (7.11) 



whereas to cut h and U we use 



h = yX" H X' H + ^(1 - y)A 3 A 3 + tA^A 3 + ^ ^ A 3 A^, (7.12) 
7 ^7 



2 



2/(1 - Z/)^h - rnm ~> 
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where 

2 

K b H = k H - ^k 3 , 7 = -sis - S23. (7.13) 

7 

All other cuts can be obtained from permutations of the ones given. It should be noted that 
the Higgs can attach itself to any gluon line, so there are quite a few more cuts than there 
are for a four-gluon amplitude. Fortunately many coefficients for our particular helicity 
arrangement are related through cyclic permutations of the external gluon momenta. 



B. Six-dimensional coefficients 

We explicitly calculate the box coefficients for Fig. [TUJ The quadruple cut is given by, 

^1234 = ^(-2)% 2 /l4 0) (-/?,l + ,2 + ,/f)4° } (-^,3 + ,if,/?) 

■<(-*!)-, l x , 2i, (la) 6 >[(-fi)*, li, 2i, (Za) 4 ] 



states 

1 



(h + k 3 )\s 3 12 

x {dlm 2 H + s 12 (h + hf) ((-^)6|/i|3 1 )[(-/ 3 ) fe |/i|3i]((/i)a|fe(-/3)l(/i)d] 

- {h + hf{s l2 m 2 H + /2(/ 1 + A ;3) 2 )((-/3) 6 |(/ 1 ),]((/ 1 )a|(-/3) b ](3l|(-/3)/ 1 |3 i ] 
+ ^12((/l + ^) 2 ^ + ^ 12 )(3 1 |(/ 1 ) (i ]((/ 1 ) a |3 i ]((-/3)^ 3 / 1 |(-/3) b ] 
-<4<H3)6|^llH8)t]<3l|(-/ 8 )/l|3i]<(il) a |fe(-/8)|(il)4]}. (7-14) 



After averaging over solutions, this becomes 



c 1 

^1234 



(12) (2 3) (31) 



mms 12 s 23 m 2 H - ^s 12 s 23 m% 



(7.15) 



[21 

We can then identify our scalar box coefficient (we neglect C\ since its integral vanishes) 
for this cut and use permutations to find it for the other quadruple cuts, 

r [o],ar> _ i ™h 

°4;123H - ~- 7Ts i.m ,v -512S23, 



2 (12)(23)(31) 
2 (12)(23)(31)' 



°4;1H23 ~~ _ o /i o\ /r. o\ /o i \ S 23'5l 3 , 



Ctim ~ 4(12)(™3)(31> Sl2Sl3 ' (T ' 16) 
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where the subscript indicates the cyclic ordering of the external particles. Our remaining 
coefficients are 



r [0]fiD 
L/ 3;123fl";12 



,6D 
/ 3;123ff;23 = 

,6D 
/ 3;12i?3;13 = 

,6D 
/ 3;12i?3;12 = 
,6D 
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r \2]fiD 
^3;1H23;1H 
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°2;12-ff3;13 " 



|4 0) (^13 + S 23), 

|4%i2 + s 13 ), 
§4%i2 + s 23 ), 

i4 0) (^3 + ^23), 

§4%i 2 + s 13 ), 
i4 0) (^ 2 + ^3), 

m H 

3 Sizmjj ' 
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^[2] 

Ly 3;123H;23 " 
^[2] 

°3;12-ff3;13 " 

°3;12-ff3;12 " 
^[2] 

°3;1H23;23 " 

r m 

°3;1H23;13 " 

^3;12H3;2H 

°2;123.ff;12 ~ 
^[21,60 
°2;lif23;23 " 



(0) 

S13+S23 



-A 



-A, 

. S12+S23 
3 

.jC) S13+S23 

3 

.ydW 312+513 

3 mi. 



(0) 

S12+S23 



-A 



-- -AAfsiz^ 

m H 

8^(0)^13^ 

8y4 (0)«i2£W 
3 S23™-Jj ' 



(7.17) 



where all other coefficients, including C2 ; i23# ; i23, vanish, and 



.4 



(0) 



H 



(12) (2 3) (31) 



(7.18) 



C. State-sum reduction 




2 1 

FIG. 11: The box diagram for the cut ordering "123H" with all internal gluons replaced by 

scalars. 



To reduce our coefficients to those of the FDH scheme, we again subtract twice the 
contribution of cuts where all internal gluons have been replaced by scalars (see Fig. flTT) . 
However, we need an effective Higgs two-scalar vertex. Using the same argument as in 
section IV D\ we can determine this vertex from the Higgs two-gluon factor by isolating the 
contribution due to internal gluons polarized along the sixth dimension while keeping the 
loop momenta in five dimensions. This effectively picks out the term with a metric tensor 
containing Lorentz indices of the internal gluons. This makes sense since it is the contractions 
of these metric tensors that gives rise to the dimensional dependence in a gluon loop when 
using the Feynman rules. The Higgs two-scalar vertex factor, as determined from eq. (17.21) . 
is then 

-2ifci • k 2 = -is 12 . (7.19) 
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We use eq. (I7.4p for the Higgs gluon two-scalar vertex factor, 

—iV2(ki - k 2 y, 



(7.20) 



where k\ and /c 2 are the momenta of the scalars. It should be noted that these are not 
the same factors that arise when taking the limit m t — > oo in a loop diagram for H — >■ ss 
mediated by a top-quark loop. 

A^\h, I s , 2 s ) is then simply given by eq. ( 17.191) . while for the Higgs gluon two-scalar tree 
amplitude we have 



4 Q) (H,ii,r,y 



—i 
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7U 
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s 12 s 13 



(7.21) 



Using these, we can find our coefficients for (H, T 



2 + , 3 + ) with a scalar loop, 
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where is given by eq. ( 17.181) . Subtracting 2 times these coefficients from our six- 
dimensional coefficients and performing the integrals, we arrive at our result in agreement 
with ref. 48 , 
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VIII. CONCLUSION 



Four- dimensional spinor helicity has proven to be a powerful tool over the years, both 
for analytic and numerical computations. It also exposes surprising structures in scattering 
amplitudes. This motivates combining the six- dimensional spinor-helicity scheme of Cheung 
and O'Connell with .D- dimensional generalized unitarity. In this paper we presented such 
a formalism as a means for simultaneously calculating cut-constructible and rational pieces 
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of one-loop amplitudes in QCD. We illustrated this by computing various sample one-loop 
amplitudes in QCD and Higgs physics. 

We obtain the amplitudes by finding the coefficients of a set of basis integrals by sewing 
together on-shell six-dimensional tree amplitudes and writing our answers in terms four- 
dimensional spinors through a decomposition of their six-dimensional counterparts. We 
extract coefficients from their limiting behavior, taking advantage of the analytic properties 
of the loop integrals, following by now standard procedures. Finally, to readjust the state 
sums to their proper values, we perform a necessary reduction procedure for internal particles 
with extra spin states in six dimensions. This is accomplished by subtracting twice the 
contribution of cuts with all internal gluons replaced by real scalars. To verify our formalism, 
we confirmed our sample calculations against known results. 

Numerical implementation should also be possible and would be interesting to explore 
and compare to other available methods. It should also be possible to apply these D- 
dimensional techniques to amplitudes with masses in the loops. Our interpretation of the 
extra-dimensional pieces as mass terms makes such computations natural in our formal- 
ism, though a modification of the integral basis to include the dropped pentagon and box 
integrals, as well as tadpoles, would be necessary. 

In summary, although many new tools are available, we believe six-dimensional helic- 
ity coupled with the unitarity method to be an additional powerful method for obtaining 
phenomenologically important loop amplitudes. 
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Appendix A: (4 - 2e)-DIMENSIONAL INTEGRAL BASIS 

In this appendix we provide details behind the integral basis given in section IIHI A quick 
review of the Feynman rules shows that a general one-loop color-ordered amplitude for n 
particles in D dimensions, where all internal particles are massless, can be written as 



where the numerator function N contains all information from external momenta and po- 
larization states, as well as tensor structure from the loop momenta. We restrict external 
momenta to be in four dimensions while allowing internal momenta to be in arbitrary D. 
Using D- dimensional Passarino-Veltman reduction techniques, we can then reduce eq. (lAlj) 
to a scalar integral basis with D-dimensional rational coefficients with at most a pentagon 




(Al) 
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integral |2l|, |23|, |49 



(4tt) d / 2 



Kb 



K 4 



+ J2 C ^s( D Kk 3 + E^;X 2 P)/& 2 , (A2) 

K 3 K 2 J 

where K r refers to the set of all ordered partitions of the external momenta into r distinct 
groups. Passarino-Veltman reduction techniques were first used in ref. 49] and a useful 
outline can be found in a recent review (Hof . 

We are focused on internal loop momenta in D = 4 — 2e. With all external momenta and 
polarization vectors in four dimensions, any contractions that they have with loop momenta 
result in four- dimensional objects. The dependence of the numerators can then only 

arise through contractions of the loop momentum with itself, 

l 2 = P + ? = P-fx 2 , (A3) 

where I contains the four- dimensional components of the loop momentum and I refers to 
the (— 2e)-dimensional components. We see then that any dimensional dependence of the 
numerators arises only through dependence on fi 2 . In a renormalizable gauge theory (such 
as QCD), the maximum number of powers of loop momentum appearing in the numerator 
of an n-point tensor integral is n, so the boxes can have at most a fi A while the triangles and 
bubbles can have up to a /i 2 . Tensor pentagon integrals are lost in the reduction, so we are 
left only with the scalar pentagon. Our basis for D = 4 — 2e then becomes 23 




+ E ^Ut^VJ + E ^kAkM 
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Further decomposition of our basis is possible by writing the scalar pentagon integral in 



terms of scalar box integrals |39|, 151 
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+ 0(e) 
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where the superscript (i) is used to refer to the box obtained from the pentagon by removing 
the propagator between legs i and i — 1. The form of the factors q is unimportant because 
they are absorbed in the box integral coefficients in our basis, which becomes 
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where ^a-kJ-a-kI ^ as absorbed terms due to the pentagon integral via eq. flA5j) . 5 In 



the same way as I 5 ~ € [/x ], e [/i ] is of order e and has been pushed into the 0(e) term. 

This recovers the basis in eq. (|3.ip . As mentioned in section HTT| evaluations of the scalar 
integrals can be found in ref. [3|, |39|. The integrals with f([i 2 ) 7^ 1 evaluate to 



jA-2t 
J 3 



jA-2e 

2;s 



6 ' 



(A7) 



where J. 4 2 ' r " 2 



2 -,, \p ] has an invariant mass square s flowing through its external legs (when the 
loop momenta are massive, I^ 26 ^ 2 ] contains additional terms). These can then be used 
to express the amplitude in terms of another common basis involving scalar integrals and 
rational terms Uli 




jA~2e 
K4 1 A;K 4 
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where 



(A9) 

In either basis finding the one-loop amplitude comes down to finding the rational coefficients 
of the non-vanishing integrals. 



Appendix B: INTEGRAL COEFFICIENT EXTRACTION DETAILS 

In this appendix, we expand upon our summary of coefficient integral extraction. The 
loop-momentum parametrizations are given in section IIVI but here we provide more details 
regarding the formulas for the coefficients. We follow Badger [23| directly for coefficients of 
integrals contributing to the rational piece of the amplitude, while simply taking /x 2 — > as 
appropriate to find the coefficients of the scalar integrals. 

1. Box integral coefficients 

When we perform a unitarity cut, we make the replacements, 

n _' , 2 ^(2n)5((l-K i ) 2 ), (BI) 



5 It should be noted that eq. (|A5j) is only valid for D = 4 — 2e. In higher dimensions, the pentagon is inde- 
pendent and its coefficient must be found separately. An additional extraction procedure is then required 
to remove its contribution to the box. For the integral basis presented here, however, no such extraction 
should be performed as the box is meant to include the pentagon contribution, even for amplitudes with 
many legs. 
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in the integral basis. Only certain integrals contain cuts depending on how many we per- 
form. Unitarity tells us that the cut integrands are proportional to products of on-shell tree 
amplitudes, so we must decompose these products to match their pieces with cut integrals, 
thereby determining the integral coefficients. 
The quadruple cut is given by 

= (47f)2_e 7 (£piE AlA ^ 4 (il), (B2) 

where the sum over a denotes the sum over the two solutions to the on-shell conditions. 
Only pentagon and box integrals contain quadruple cuts, so 

(4 * )2 ~'/ (S^E^^W) 

K 5 Ka K 4 K 4 

(B3) 

where -^ ut [/(/i 2 )] is the quadruple-cut version of the (4 — 2e)-dimensional integral in eq. ( 13. 2j) . 
and a factor of | has been included to account for the averaging over solutions. In simple 
cases, the integrand is a polynomial in fi 2 . This is the case for four-point amplitudes with 
massless particles, which also do not contain a pentagon contribution. After integrating the 
cut integrals over the delta functions, we can easily identify our coefficients by matching 
orders of fi 2 in 

Cf + fcf + fcf = \ £ A 1 A 2 A,A4l). (B4) 

cr 

Even in cases where we have a nonzero C 5 , it would not have to be disentangled from C\ 
since the integrals associated with them are of order e. In more general cases, the quadruple 
cut does not fall naturally into a polynomial in fi 2 , and we must identify each coefficient. 
We can find Cf , our cut-constructible piece, by taking /z 2 — > 0, 

which effectively mimics a four-dimensional generalized unitarity procedure. To find Cf, 
following ref. [23(, we express our integral with a term containing the behavior at infinity 
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and residue terms contained in a sum over poles in Fi 

= {4n) / WF* E [Inf, 2 AA 2 A 3 A 4 ^)](/i 2 ) + £ o_ 2 

^ ' CT \ poles{i} " ^ l 

(B6) 

As mentioned in section IIV} the Inf term can be expressed as a polynomial in F with 
maximum order Fi 

2 

[M^A^A^F) = J] q/i 2 \ (B7) 

i=0 

To express the residue terms as a polynomial in /i 2 , we partial fraction and find at most a 
term of order /i 2 , 



Res^2(A 1 A 2 A 3 A 4 (l^)) _ fx 2 Res^(A 1 A 2 A 3 A 4 m)) Res ^ = ^ ( A, A 2 A 3 A 4 (l^)) 



(B8) 



v 2 - F F(v 2 - F) F 

Evidently the first term in eq. flB8j) combines with the order F term in the boundary 
polynomial to give C\ entangled with the cut pentagon, and the second term combines 
with the order /i° term to give Cf\ cjj ' could be found by combining these terms, but 



eq. (IB5I) is much simpler. The crucial point is that only the boundary term has a term of 



order /i , so we find 



C f = I E M^AAsM/^L 4 , (B9) 



where the Inf operator has been restricted to the coefficient of the F term. Along with 
eq. (1B5I) . we have the box coefficient extraction formulas in eq. (14. 8p . 

2. Triangle integral coefficients 

Recall that the triple-cut momentum solution has a free complex parameter t, 

11 = aKf + bK* + UkI^KI) + {Kl\r\Kl). (BIO) 

We can then express a generic triple cut as 

(4 ! r) 2 -<(-2 I ri) s | ^^W^A^A, 



i - 1 7 - 

poles {j } 



(Bll) 
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In going from the second line to the last line, we have performed the integral transformation 
from to t and integrated over the three delta functions (thereby applying the on-shell 
constraints which give rise to the sum over solutions, including solutions for 713 and the 
conjugate momentum solution). The Jacobian of the transformation is J t (which also absorbs 
a factor of l/2ir). As before, the Inf term contains the information from the boundary of 
the t contour integral, 

\im([lni t A 1 A 2 A 3 }(t)-A 1 (t)A 2 (t)A 3 (t))=0 J (B12) 

t— >oo 

while the second term is a sum over poles in t. 

Poles in t in the cut integrand come from propagator terms of the form 1/(1 — K) 2 . The 
only contributions from the residue terms occur where t = tj, which corresponds to putting 
another propagator on shell (I — K) 2 = (I — K) 2 — fi 2 = 0. Therefore the second term 
contains only box and pentagon coefficients, and we can ignore it in our search for triangle 
coefficients, 

= z(47r) 2 ~ e J J dtJ t J2l^tA 1 A 2 A 3 (F 1 )](t) + box/pent, terms. (B13) 

For the first term, much like before we can expand around infinity and write it as a 
polynomial in t, 

m 

[lni t A 1 A 2 A 3 }(t) = Y,tt i - (B14) 

Because of our particular parametrization of the loop momentum fIBlOl) . all integrals over 
t n for n 7^ vanish (§1, 



J dtJ t t n = 0, n ^ 0, (B15) 
allowing us to drop all terms in the expansion not of order t°. Our triple cut is then 

= i(4n) 2 ~ e J ^-y^ hVi) j dtJ t + box/pent, terms 

= i(4ir) 2 - e J AJLj^ [[lni t A 1 A 2 A 3 (ll)}(t)\ t ^ ] J dtJ t + box/pent, terms. 

(B16) 

The triple-cut triangle integral is precisely 

-(4-)^/ (£p^ 2 ) / dtJ^i^r^ifj *^ftf)f[8(!f), (B17) 
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which allows us to identify the triangle coefficients and the appropriate sign. For simple 
cases where the Inf term falls into a polynomial in /i 2 , we can match orders of \i. Averaging 
over solutions (both of 713 and the conjugate momentum solution for fixed 713), we have 

Cf + tfCf = ~Y, [Inf^^lWLo , (B18) 

where n 7 denotes the number of solutions for 7 13 , which is either 1 or 2. In more general 
cases where [lnf t AiA 2 A 3 (l^)](t) | is not a polynomial in /1 2 , we use a procedure similar to 

that of section IB II We take /i 2 — > for Cf 1 , and for Cf we use a polynomial expansion in 
tr- 

1 

[Inf M2 [Inf^ 1 A 2 A 3 (/l)](t)](/. 2 )|^ = $>/i 2 \ (B19) 

i=Q 

and restrict to the coefficient of the fj, 2 term, recovering eq. (14.131) . 

^ a 

C f = [Inf, 2 [Inf^ 1 A 2 A 3 (/n](t)](^ 2 )| M2 ^ • (B20) 



3. Bubble integral coefficients 



The double-cut momentum solution (I4.14p has two free parameters t and y, so the double 
cut is given by 
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(B21) 



where we have integrated over the delta functions and performed an integral transformation 
from Z M to t and y, similar to our treatment of the triple cut. The last term of the final 
expression has two additional propagators on shell, and its numerator has no dependence 
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on t or y. It therefore corresponds to box and pentagon terms only. We might then believe 
that the second and third terms correspond to triangle coefficients only and grab the first 
term for the bubbles, but the fact that there is dependence on the free parameters t and y 
in the numerators is significant. 

These single residue terms have a third propagator on shell, which indeed allows us 
to correspond them to a triple cut. We can, without loss of generality, fix y to satisfy the 
constraint and use a triple cut to find the residue terms rather than computing them directly. 
However, this parametrization of the loop momentum may differ from that of section IB 21 
depending on our choice of X . With the new parametrization, integrals over positive powers 
of t that conveniently vanished in section lB~2l will not necessarily vanish here, 



J dt,J' t t n ^ 0, (B22) 



where J' t is the Jacobian corresponding our new loop momentum parametrization. This 
effectively means that there are powers of the loop momentum (containing more than just 
the extra-dimensional pieces /x) in the numerators, and our triple cut produces tensor triangle 
integrals. These tensor integrals are not in our integral basis and must be reduced. Passarino- 
Veltman reduction shows us that there are both scalar triangle and scalar bubble integrals 
within these tensor integrals, 



dtJ[t n = C^T + C 2 IT , (B23) 

so we must extract the bubble integral coefficients. 

Imposing a third on-shell condition in our momentum parametrization, eq. (14.141) gives 
for y, 



y± = Wf + 4 B ^ (B24) 



where 



B 2 = S 1 ( X -\$^\K\-), 

B 1 = jt{Ki~\#: 3 \Kt) ~ S 1 t( X -\^ 3 \x') + S 1 ( X -\^ 3 \Ki-), 

B = ^t 2 (K b 1 -\^ 3 \x-)-^(x'\^\K b 1 -)+itS 3 + tS 1 ( X -\^ 3 \x-)- (B25) 

After dropping box and pentagon terms, the triple cut is given by 

/ dtJ'AlnUA^A^m, (B26) 



where the sum is over the solutions in eq. flB24j) . Once again we can expand around t — > 00 



[lnf t A 1 A 2 A 3 ](t) = Y,M\ (B27) 



8=0 



to perform our integrals. The t° term gives us the double-cut scalar triangle, while positive 
powers of t return, among other pieces, cut scalar bubble coefficients. We evaluate the 
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integrals over positive powers of t and retain only the contributing bubble integral to our 
particular double cut, 

(4tt) 2 j dtJ[t j = T,/ 2 4;cut 

= ~ 1? J * { h Cj 'a^-i^ 1 J /s ' (B28) 



where 
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For QCD we can have terms up to order t 3 , so our relevant terms are 

SiOrlJfsltfJ-) 
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(B30) 

recovering eq. f)4.20p . We then have the contribution of the single residue terms to the bubble 
coefficient. Defining T = 0, we have 

-ij2[lnhA 1 A 2 A 3 ](t)\^ Tj , (B31) 

Uy 

after averaging over solutions. Single residue terms arise from every possible third leg that 
can go on-shell, i.e. every triple cut that shares two of its cuts with our two-particle cut, so 
we must find a term for every such possible configuration. 

Lastly we need the bubble coefficients due to the first term in eq. (1B21j) . 

-(4^| -J^ J dtdyJ t , y J2[m^yAiA 2 (imy)}(t)- (B32) 
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Once again we expand around infinity to evaluate these integrals. Because of our choice of 
momentum parametrization, integrals of positive order in t vanish. This is not the case with 
y, so we must keep terms of the form t°y m . These evaluate to 



Y m = (2tt) 2 J dtdyJ t , y y T 
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in which we identify the double-cut scalar bubble integral, 

i{^f- e J J dtdy.J t , y = {-i){Anf-%-2nif j ^^f[ 5 ^- (B36) 

Adjusting a factor of i, we can easily identify the bubble coefficient from the first term of 
eq. flB21~]) as 

-i [Inf t [Inf „ A 2 ] (j/)] (*) |*->o l2 r->y m • (B37) 

We then have for our total bubble coefficient, in the case where our calculation falls into a 
polynomial in /i 2 , 

C? +^Cf = -i[InU[lnf y A 1 A 2 ] (y)] {t)\ t ^ y m^ Ym - ^ ^prf^^]®^, 

Ctri 

(B38) 

where Cm denotes a sum over all possible triangles attainable from cutting one more leg of 
our two-particle cut. Most generally, we have 

C l ° ] =-i\M t \M y A 1 A 2 ](y)](t)\ li 

->0,£> ->Tj ) 
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in agreement with eq. (I4.18p . 
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